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The properties of the momentum space quark propagator in Landau gauge are examined for the 
overlap quark action in quenched lattice QCD. Numerical calculations are done on three lattices with 
different lattice spacings and similar physical volumes to explore the approach of the quark propa- 
gator toward the continuum limit. We have calculated the nonperturbative momentum-dependent 
wave function renormalization function Z{Q 2 ;p) and the nonperturbative mass function M(p) for a 
variety of bare quark masses and perform an extrapolation to the chiral limit. We find the behavior 
of Z(( 2 ;p) and M(p) are in reasonable agreement between the two finer lattices in the chiral limit, 
however the data suggest that an even finer lattice is desirable. The large momentum behavior is 
examined to determine the quark condensate. 
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>. ' I. INTRODUCTION 

00 ' 

The quark propagator is one of the fundamental quantities in QCD. By studying the momentum-dependent quark 
mass function in the infrared region we can gain valuable insight into the mechanism of dynamical chiral symmetry 
breaking and the associated dynamical generation of mass. There have been several studies of the momentum space 
£f) . quark propagator Q, 0, 0, 0, 0, IE IE El m Landau gauge using different fermion actions. Here we focus on the 
overlap fermion action and extend previous work § to three lattices with different lattice spacing a at fixed physical 
volume. This allows us to study the approach of the Landau gauge quark propagator to the continuum limit. The 
study of the overlap quark propagator in the Gribov copy free Laplacian gauge is underway and will be reported 
elsewhere. 
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II. QUARK PROPAGATOR ON THE LATTICE 

In a covariant gauge in the continuum, the renormalized Euclidean space quark propagator has the form 



<D 
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X 

■- : s(e-p)= 1 - z(c2;p2) a) 

V ^ ,P) ij>A(( 2 ;p 2 )+B(( 2 ;p 2 ) z^ + M(p 2 )' V > 

where £ is the renormalization point. The renormalization point boundary conditions are chosen to be 

Z(( 2 ;( 2 )^l, M(C 2 )^m(C 2 ). (2) 

where m(£ 2 ) is the renormalized quark mass an ultraviolet renormalization point. The functions A(£ 2 ;p 2 ) and 
B(( 2 ;p 2 ), or alternatively Z(( 2 ;p 2 ) and M(p 2 ), contain all of the nonperturbative information of the quark propagator. 
Note that M(p 2 ) is renormalization point independent, i.e., since S(C, 2 ;p) is multiplicatively renormalizable all of the 
renormalization-point dependence is carried by Z((^ 2 ;p 2 ). For sufficiently large momenta the effects of dynamical chiral 
symmetry breaking become negligible, i.e., for large p 2 and we have M(p 2 ) — > m(£) up to logarithmic corrections, 
where m(C) is the perturbative running mass. 

When all interactions for the quarks are turned off, i.e., when the gluon field vanishes (or the links are set to one), 
the quark propagator has its tree-level form 

Si0) ^=wrr^> ( 3 ) 

vp + m 
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where m° is the bare quark mass. When the interactions with the gluon field are turned on we have 

S<°>(p) -> S baro (a;p) = Z 2 (C 2 ;a)5(C 2 ;p) , (4) 

where a is the regularization parameter - in this case, the lattice spacing - and Z2(C, 2 ]o) is the quark wave- function 
renormalization constant chosen so as to ensure Z(£ 2 ;p 2 ) = 1. For simplicity of notation we suppress the a-dependence 
of the bare quantities. 

On the lattice we expect the bare quark propagators, in momentum space, to have a similar form as in the continuum, 
except that the 0(4) invariance is replaced by a 4-dimensional hypercubic symmetry on an isotropic lattice. Hence, 
the inverse lattice bare quark propagator takes the general form 



(s^r^^i^c^^j +b( P ). 



(5) 



We use periodic boundary conditions in the spatial directions and anti-periodic in the time direction. The discrete 
momentum values for a lattice of size Nf x N t , with n.$ = 1, .., Ni and nt — 1, .., N t , are 

271 ( N A A 27T (m 1 N A (fi\ 

Pi= N^{ ni -T)> ^ Pt ^Wa{ Nt ~2~Y)- (6) 

The overlap fermion formalism |lOtlllj| realizes an exact chiral symmetry on the lattice and is automatically O(a) 
improved. The massive overlap operator can be written as [12j 

D(ri = ^[l + [i+(l-fi) l5 e(H w )} , (7) 

where H w (x, y) = y^ ) D w {x 1 y) is the Hermitian Wilson-Dirac operator, the mean-field improved Wilson-Dirac operator 
can be written as 



D w (x, y) = [{-m w a) + 4r] S XjV - i ^ {(r - y^)U fl (x)S VtX+ / 1 + (r + y^Ufa - afi)5y tX -p,} 

(8) 
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The negative Wilson mass (—m w a) is then related to k by 

u 



2(—m w a) + (1/k c ) 



(9) 



and mean- field improvement allows the use of the tree-level value k c = l/(8r). The Wilson parameter is typically 
chosen to be r = 1 and we will also use r = 1 here in our numerical simulations. The dimensionless quark mass 
parameter is 

(10) 

2m w 

The overlap quark propagator is given by the equation 

5 bare (m°) = D~ x (n) , (11) 

where 

D-\n) = ^D-\fj) and D~ x (n) = ^-^ — [-D _l (Ai) — l] . (12) 

At tree-level, the inverse bare lattice quark propagator becomes the tree-level version of (Eq. |5J) 



(13) 
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We calculate S^°\p) directly by setting the links to unity in the coordinate space quark propagator and taking its 
Fourier transform. It is then possible to identify the appropriate kinematic lattice momentum directly from the 
definition 

Q* = CW(p) (14) 

The form of q^ip^) is shown and its analytic form given in Ref. Having identified the appropriate kinematical 
lattice momentum q, we can now define the bare lattice propagator as 

This ensures that the free lattice propagator is identical to the free continuum propagator. Due to asymptotic freedom 
the lattice propagator will also approach the continuum form at large momentum. In the gauge sector, this analysis 
approach dramatically improves the giuon propagator [l3lll4). 
The two Lorentz invariants can now be |25j 

Z- 1 (p) = j ^TtUS- 1 (p)} (16) 

M(p) = ^Tr{5- 1 (p)}. (17) 

While Z is directly dependent on our choice of momentum, q, the mass function M is indirectly dependent on this 
choice. In the case of staggered quarks it has been seen that the kinematic momentum derived from tree-level analysis 
of the action is a good choice of momentum for the mass function 0,0. This is an empirical result. The tree- level 
behavior of the Overlap quark propagator is rather different, however, and a different approach may be needed. We 
investigate this issue by analyzing the scaling behavior of the propagator over three values of the lattice spacing at 
constant physical volume. 



III. NUMERICAL RESULTS 



We present results from three lattice ensembles, each with a different lattice spacing, a, but having the same physical 
volume. Lattice parameters are summarized in Table [I] The gauge configurations were created using a tadpole 
improved plaquette plus rectangle (Liischer-Weisz 0) gauge action. Each ensemble consists of 50 configurations. 
The lattice spacing was determined by the static quark potential using the string tension \fa = 440 MeV [16( . 



TABLE I: Lattice parameters. 



Action 


Volume 


AT-Therm 




P 


a (fm) 




Physical Volume (fm 4 ) 


Improved 


16 3 x 32 


5000 


500 


4.80 


0.093 


0.89650 


1.5 3 x 3.00 


Improved 


12 3 x 24 


5000 


500 


4.60 


0.123 


0.88888 


1.5 3 x 3.00 


Improved 


8 3 x 16 


5000 


500 


4.286 


0.190 


0.87209 


1.5 3 x 3.00 



The gauge field configurations were gauge- fixed to the 0(a 2 ) improved Landau gauge 18]. Our calculation begins 
with the evaluation of the inverse of the Dirac operator in Eq. Q. We approximate the matrix sign function e(H w ) 
by a 14th-order Zolotarev approximation [l^. The coordinate space propagator, Eq. I|ll|) . is calculated for each 
configuration. A discrete Fourier transform is then applied to the each of the coordinate space propagator, and the 
momentum-space bare quark propagator, 5' bare (p), is finally obtained from the ensemble average. 

In the Wilson action we use n — 0.19163 for the regulator mass. We calculate the Overlap quark propagator for 
ten quark masses on each ensemble by using a shifted Conjugate Gradient solver. The quark mass parameter \i was 
adjusted to make the tree level bare quark mass in physical units, the same on three lattices. For example, we choose 
fi = 0.018, 0.021, 0.024, 0.030, 0.036, 0.045, 0.060, 0.075, 0.090, and 0.105 on ensemble 1, i.e., the 16 3 x 32 lattice 
with a = 0.093 fm. This corresponds to bare masses in physical units of m° = 2p,m w — 127, 148, 169, 211, 254, 317, 
423, 529, 634, and 740 MeV respectively. 

Results from ensemble 2 were presented in Ref. § , and some results from ensemble 3 were also reported in Ref. [i| . 
Here we will compare the quark propagators on each ensemble to examine its behavior as the lattice spacing vanishes. 
First we present some results from ensemble 1, the finest lattice of the three. All data has been cylinder cut [l3| . 
Statistical uncertainties are estimated via a second-order, single-elimination jackknife. 
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FIG. 1: The functions M{p) and Z {K) (p) = Z{C, 2 ;p) renormali zed at ( = 3.44 GeV for all ten quark masses. Data are plotted 
versus the discrete momentum values defined in Eq. (JSJ, p = x/SPh; over the interval [0,5] GeV. The data correspond to bare 
quark masses (from bottom to top) fi = 0.018, 0.021, 0.024, 0.030, 0.036, 0.045, 0.060, 0.075, 0.090, and 0.105, which in physical 
units corresponds to m° = 2 f j,m w ~ 127, 148, 169, 211, 254, 317, 423, 529, 634, and 740 MeV respectively. 



In Fig. ^ w e show the results for all ten masses for both the mass and wave function renormalization functions, 
M(p) and Z^ R \p) = Z{C, 2 ;p) respectively, as a function of the discrete lattice momentum p. Z^ K \p) is renormalized 
at £ = 3.44 GeV. We see that both M(p) and Z^{p) are reasonably well-behaved up to 5 GeV. In the plots of M(p) 
the data is ordered as one would expect by the values for fj,, i.e., the larger the bare quark mass m° the higher is the 
M(p) curve. At small bare masses M(p) falls off more rapidly with momentum, which is understood from the fact 
that a larger proportion of the infrared mass is due to dynamical chiral symmetry breaking at small quark masses. 
In the non-relativistic limit the mass function would be a constant. 

Z^ RS> (p) on the other hand is infrared suppressed. The smaller the quark mass, the more pronounced the di p at low 
momenta. This behavior is qualitatively consistent with what is seen in Dyson-Schwinger based QCD models|20ll2l|. 
It is likely that some of the suppression, however, is due to the finite volume 6]. In Fig. [5] we repeat these plots but 
now using the kinematical lattice momentum q. This only alters the large momentum behavior of the propagator. 

We perform an extrapolation to zero quark mass by a linear fit to the data. At sufficiently large momenta the mass 
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FIG. 2: The functions M(p) and Z^ K \p) = Z((^ 2 ;p) for renormalization point £ =5.31 GeV for all ten quark masses. Data are 
shown versus the discrete momentum values defined in Eq. 114H . q = \/Yl 9ui over the interval [0,12] GeV. The data in both 
parts of the figure correspond from bottom to top to increasing quark masses. The values of the bare quark masses are in the 
caption of Fig. Q 



function will be proportional to the bare quark mass, in which case the linear extrapolation is appropriate. Non- 
linear behavior is to be expected in the infrared, but this simple ansatz describes the present data adequately. In the 
ultraviolet, the renormalized Z should - and does - approach its perturbative value of 1. This is mass independent. We 
investigated the effect of including the smallest quark masses in the chiral extrapolation and found that eliminating the 
lightest two made little difference to the extrapolated result. The resulting estimate of the chiral limit is shown in Fig. 13 
These are shown against both p and q, renormalized as before. We see that both M(p) and Z^ K \p) deviate strongly 
from their tree- level behavior. In particular, as in earlier studies of the Landau gauge quark propagator00,IE@l) we 
find a clear signal of dynamical mass generation and a significant infrared suppression of the Z^ R \p) function. At the 
most infrared point - the lowest non-zero momentum available on this lattice - the dynamically generated mass has 
the value Mir = 307(6) MeV and the momentum-dependent wave function renormalization function has the value 
Zm. = 0.55(2). These values are very similar to the results found in previous studies @i IH IS II] and are also similar 
to typical values in QCD-inspired Dyson-Schwinger equation models [2D, EI- The results of Ref. suggests that at 
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FIG. 3: Linearly extrapolated estimates of M(p) and Z^ K \p) = Z(£ 2 ;p) in the chiral limit. Here the renormalization point are 
£ =3.44 GeV in the p scale and £ =5.31 GeV in the q scale. At the smallest accessible momentum Mir = 307(6) MeV and 
Zir = 0.55(2). 



least some of the infrared suppression of Z^ R \p) is due to finite volume effects. 

Now we present the results on three lattices for comparison. These lattices have approximately the same physical 
volume, but each has a different lattice spacing. Thus we can study the Overlap propagator's scaling properties. We 
present the results for the chiral limit. The mass function, M(p) for the three lattices in the chiral limit is plotted in 
Fig. using both p and q. We see that if the mass function M(p) is plotted against the standard lattice momentum 
p, the agreement of the results among the three lattices is better than the case in which the mass function M(p) is 
plotted against kinematical lattice momentum q. 

The results for the renormalization function Z^(p) = Z(( 2 ;p) of the three lattices is plotted in Fig. [5] Here 
the renormalization point are chosen to be £ = 3.44 GeV in p scale and £ = 5.31 GeV in q scale. Contrary to the 
case of mass function M(p), but as predicted by the tree-level analysis, the agreement between the results on three 
lattices is better if Z^(p) is plotted against the kinematical lattice momentum q. There are also the relatively small 
discrepancies in Z^\p) versus q in the infrared region on three lattices, it seems that in the continuum limit, the 
dip in the renormalization function Z^ljp) will be narrow but the depth of the dip will be unchanged. It suggests 
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FIG. 4: The mass function M(p) from a linear extrapolation to the chiral limit is shown for our three lattices. In the upper 
part of the figure M(p) is plotted against the discrete lattice momentum p, whereas in the lower part it is plotted against the 
kinematical momentum q. The results again suggest that we most rapidly approach the continuum limit by plotting M(p) 
against p. 



that an even finer lattice will be needed to confirm the continuum limit of (p) in the infrared. It is possible that 
the linear chiral extrapolation is unreliable for Z^ R ' (p) in this regime or it could be that dynamical chiral symmetry 
breaking is coupling hypercubic lattice artifacts to finite volume effects. This warrants further investigation with finer 
and larger lattices. 

Thus we have resolved one of the key questions raised in the studies of Ref. @ . We see that the continuum limit 
appears to be approached most rapidly when Z^{p) is plotted against q and M(p) is plotted against p. The better 
scaling of Z^ K '{p) as a function of q is natural and predicted by the tree-level analysis. That M{p) is better against 
p is purely observation. 

Another way of studying scaling is by making comparisons with known continuum results. In the chiral limit, in 
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FIG. 5: The momentum-dependent wave function renormalization function Z^(p) = Z(( 2 ;p) for renormalization point 
£ = 3.44 GeV in the p-scale and £ = 5.31 GeV in the g-scale from a linear extrapolation to the chiral limit. In the upper 
part of the figure (p) is plotted against the discrete lattice momentum p whereas in the lower part it is plotted against the 
kinematical momentum q. The results again suggest that we most rapidly approach the continuum limit by plotting Z^'(p) 
against q. 



the continuum, the asymptotic quark mass function has the form 

r -^_4^ m w^'i'- (18) 

3 [ln(/i 2 /A| CD )] d M p 2 

(see Ref. [2(j, Eq. (6.15)) where the anomalous dimension of the quark mass is cLm = 12/(33 — 2JV/) for Nf quark 
flavors (Nf = in the present case). The dependence of M(p 2 ) on the renormalization point /.i is canceled by the 
dependence of the condensate, maintaining the renormalization point invariance of the mass function. We fit this 
form to the lattice data obtained by both linear and quadratic chiral extrapolation. A quadratic extrapolation was 
used on the AsqTad data [(| so it is useful for comparing with those results. Some results are presented in Table ITU 
The difference between the quadratic and linear extrapolations is no great surprise as our quark masses are rather 
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TABLE II: Extracted values of the quark condensate 



p 


extrapolation 


p fit 


region (GeV) 


-{^V> 1/3 (MeV) 


q fit 


region (GeV) 


-(^V) 1/3 (MeV) 


4.60 


linear 


3.6 - 


4.5 


337(39) 


4.3 - 


8.6 


621(49) 


4.60 


quadratic 


3.6 - 


4.5 


292(56) 


4.3 - 


8.6 


575(72) 


4.80 


linear 


3.6 - 


5.3 


327(22) 


5.5 - 


11.4 


499(34) 


4.80 


quadratic 


3.6 - 


5.3 


259(36) 


5.5 - 


11.4 


395(54) 



heavy. This is a constraint of the volume. The relevant point is that there is good agreement between the two lattices 
when p is the momentum, but not when q is the momentum. 

IV. SUMMARY AND OUTLOOK 

The momentum space quark propagator has been calculated in Landau gauge on three lattices with different lattice 
spacing a but very similar physical volumes in order to explore the approach to the continuum limit. We calculated 
the nonperturbative momentum-dependent wave function renormalization Z(£ 2 ;p) and the nonperturbative mass 
function M(p) for a variety of bare quark masses. We also explored the quark propagator in the chiral limit. As 
previously anticipated [8|, the continuum limit for Z(C, 2 ;p) is approached most rapidly when it is plotted against the 
kinematical lattice momentum q, whereas for the quark mass function, M{p), we have found that using the discrete 
lattice momentum p provides the most rapid approach to the continuum limit. 

Future work should test our conclusions and further explore the continuum limit with one or more additional finer 
lattice spacings. In addition, it will be necessary to use both finer and larger volume lattices, in particular to study 
the infrared behavior of Z^ 2 , q). One can also use other kernels in the overlap fermion formalism, e.g., using a fat-link 
irrelevant clover (FLIC) action (2^ as the overlap kernel [23l |24| in order to further establish the robustness of our 
conclusions and to provide more accurate data. These studies are currently underway and results will be reported 
elsewhere. 
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